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Abstract 

This work explores the relationship between optimal control theory 
and adiabatic passage techniques in quantum systems. The study is based 
on a geometric analysis of the Hamiltonian dynamics constructed from 
the Pontryagin Maximum Principle. In a three-level quantum system, 
we show that the Stimulated Raman Adiabatic Passage technique can 
be associated to a peculiar Hamiltonian singularity. One deduces that 
the adiabatic pulse is solution of the optimal control problem only for a 
specific cost functional. This analysis is extended to the case of a four-level 
quantum system. 

1 Introduction 

Adiabatic passage techniques are an efficient approacli to design a control pulse 
wfiicfi drives a quantum system from an initial state to a specific target state (for 
recent overviews, see [TJ [2] and references tlierein) . Sucfi processes liave found 
a large variety of applications ranging from Nuclear Magnetic Resonance [3] 
and quantum information to atomic and molecular excitations [H [5] . Adiabatic 
metliods are usually achieved by using a series of intense pulses which can be 
frequency chirped, the frequencies and the chirped being chosen with respect 
to the structure of the energy levels of the quantum system. The modification 
of the shape of the pulse envelope and the chirping rate must be sufficiently 
slow so as to fulfill adiabatic conditions. In this paper, we will consider one 
of the most well-known adiabatic processes, the Stimulated Raman Adiabatic 
Passage (STIRAP) excitation [6l[7], which involves a counterintuitive sequence 
of two pulses in a three-level quantum system. This approach, which allows a 
complete population transfer between two energy levels, is also robust in the 
sense that it is not sensitive to small variations of the field parameters. The 
adiabatic methods have however some major drawbacks, which can be roughly 
summarized in a high pulse energy and a long control duration. Such constraints 
can be problematic, e.g., if other concurrent physical or chemical processes with 
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the same time scale occur during the control. This statement can be expressed 
by saying that adiabatic fields are not able to reach the physical limits of the 
best possible performance of the system in terms of duration or energy of the 
pulse. These limits can be established by using optimal control theory [HlElllllS]) 
which has been introduced in the eighties in the quantum control community. 
An optimal control problem can be solved either from geometric and analytic 
tools [ini [TTJ [ini [ini [TSI IU] tor systems of low dimensions or from numerical 
algorithms in the case of large dimensional quantum systems [SI [51 [T7[ [TSj . 
The two approaches are based on the application of the Pontryagin Maximum 
Principle (PMP) [191 [2T], which transforms the control problem into an 
Haniiltonian one. In this framework, the computation of the optimal solution 
is replaced by the determination of an Hamiltonian trajectory satisfying given 
boundary conditions. A general problem with optimal solutions is their lack of 
robustness against variations of the system parameters and of the control field. 
Note that a method based on the simultaneous control of different systems can 
be used to improve this latter property (see e.g. [22]), but this purely numerical 
approach does not give any insight into the robustness nature of the optimal 
solutions. 

In the present paper, we study the relationship between adiabatic and op- 
timal pulses in the case of a three-level quantum system where the STIRAP 
solution can be used. This question has been analyzed in different works over 
the past few years but not directly from a geometric optimal control perspec- 
tive. This work complements therefore several earlier numerical [531 US [551 1211 
[ini [571 [551 [Ml [Sni 131] and analytical studies [Ml [10] on this relation. To sum- 
marize some previous results, we mention that the authors of [27] argue that 
the STIRAP technique cannot be expressed as an optimal control solution. The 
counterintuitive sequence of STIRAP was recovered numerically in |26| by con- 
sidering as cost functional the population of the intermediate level. In Ref. [TD] , 
a tracking technique aiming at minimizing the population into the intermediate 
state was proposed to recover the adiabatic solution, but no insight into the 
relation between the adiabatic pulse and the optimal one was given. In [32], 
the Hamilton- Jacobi Bellman theorem was used to solve the optimal control 
problem. 

In this work, we first consider the energy minimization control problem. Us- 
ing the Hamiltonian structure derived from the PMP, we show that the STIRAP 
solution can be associated to a peculiar singularity of this system. However, in 
this case, the adiabatic pulse can be viewed as a singular limit solution of the 
Hamiltonian system which cannot be approached smoothly by an optimal con- 
trol field. In a second step, we show that the STIRAP pulse is solution of an 
optimal control problem corresponding to a specific cost functional, that we call 
STIRAP cost. Using this geometric analysis, we extend this approach to the 
case of a four-level quantum system [3 [331 • 

The paper is organized as follows. We present the system in Sec. 12. II and we 
recall some basics about the STIRAP solution and the application of the PMP 
in Sec. 12.21 Section [5?31 is devoted to the geometric analysis of the Hamiltonian 
dynamics for the energy minimization problem. It is shown in Sec. 12.41 that 
the adiabatic pulse is solution of an optimal control problem for a specific cost 
functional. In Sec. [31 we extend this study to the case of a four-level quantum 
system. A summary is presented in Sec. 31 The appendix [^ collects some 
technical and mathematical results. 
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2 Optimal control of a three-level quantum sys- 
tem 



2.1 The model system 

We consider a three-level quantum A— system whose dynamics is governed by the 
Schrodinger equation. The system is described by a pure state \tp{t)) belonging 
to a three-dimensional Hilbert space 7i spanned by the basis {|1), |2}, |3)}. The 
dynamics of the system is controlled by the pump and the Stokes pulses which 
couple respectively the states |1) and |2) and the states |2) and |3). Note that 
there is no direct coupling between the levels |1) and |3). The time evolution of 
\tjj{t)) is given by 

^^^\^it))^Hit)\m), (1) 

where, in the interaction representation and the rotating-wave approximation, 
the Hamiltonian H{t) can be written as 

/ ui{t) \ 
H{t) = u,{t) -ik -U2{t) , (2) 
V -U2(t) / 

with ui(t) and U2(t) representing the Rabi frequencies of the pump and Stokes 
pulses, respectively. The two pulses are assumed to be on-resonance with the 
corresponding frequency transitions. Equation ([T]) is written in units such that 
h= 1. The parameter k describes the relaxation rate of the second level, which 
is the only state interacting with the environment. Its role will be made clearer 
in the following. Note that the Hamiltonian ^ has already been considered 
in other studies about the use of adiabatic passage techniques in presence of 
dissipation [321 [551 137j . We denote by ci, C2 and C3 the complex coefficients 
of the state li'it)), and we introduce the real coefhcients Xi {i G {1, 2, • • • , 6}) 
defined by: 

ci = xi + 1x4, C2 ^ X5- ix2, C3 ^ X3 + ixe- (3) 

Using Eq. ([T]). it is then straightforward to see that only the coordinates xi, X2 
and X3 are coupled to each other. This leads to the following differential system 

Xi \ I -U1X2 \ 

±2 = UiXi - kX2 + U2X3 , (4) 

X3 / \ -U2X2 J 

which reads in a more compact form as: 

X = Fo{x) + uiFi{x) + U2F2{x), (5) 

with X = {xi,X2,X3), Fo{x) = (0, -fca;2,0), Fi{x) = (-.T2,a;i,0) and F2{x) = 
(0,3:3, -X2). 

We aim at transferring the system from the state |1) to the state |3) without 
losing population through the dissipation of the state |2). In the new coordi- 
nates, this corresponds to the passage for xi = 1 to X3 = 1. This transfer can 
be realized by a standard STIRAP strategy in which the pulses are applied in 
a counterintuitive order, i.e. the Stokes pulse precedes the pump pulse. We 
propose to revisit in the next section this control problem by using the PMP. 
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2.2 The Pontryagin Maximum Principle 

We first analyze the optimal control of this three-level system with the constraint 
of minimizing the energy of the laser fields. The total duration T is fixed 
and there is no restriction on the amplitudes of ui and U2. For this energy 
minimization problem, the cost functional C can be written as: 

r[ul{t)+ul{t)]dt, (6) 



where T is the fixed control duration. The PMP states that there exist an 
adjoint state p — {pi,P2,P3) and a negative constant pq not simultaneously equal 
to zero such that the optimal trajectory is solution of the following Hamiltonian 
problem: 

H^p-{Fo+ ELi P^) +Po{ul{t) + ul{t)) 
x=^{x,p,v) 

H{x,p,v) — max„iJ(x, p, m), 

where '•' stands for the scalar product between two vectors. This optimal control 
problem with no dissipation, i.e. k = 0, has been already solved in Ref. |10) 
where it was shown that the optimal solution is associated to an intuitive order 
of the Stokes and pump pulses. Here, we consider the same analysis by adding 
a dissipative term on the intermediate state which forces the control to not 
populate this level, mimicking thus the adiabatic trajectory. Since there is no 
constraint on the control fields, these latter can be computed explicitly from the 
maximization condition which leads to: 

— - — - (8) 
dui ' du2 

In the regular case where pq can be normalized to —1/2, one deduces that the 
optimal controls are given by ui — p-Fi and U2 — p-F2. Plugging the expressions 
of ui and U2 into the first equation of it is straightforward to check that 
the Hamiltonian H reads: 

H=p-Fo + ^{p-Fi)^ + ^ip-F2)^ (9) 

The corresponding Hamiltonian dynamics encodes all the information about the 
optimal trajectories of the control problem. A global overview of these solutions 
can be obtained through a geometric analysis of this dynamics, which is made 
clearer by the introduction of the following spherical coordinates: 

: r sin 9 cos (f) 

= rcos6l (10) 
: r sin 9 sin 0. 

The dynamical system takes then the form: 

= —krcos^O 

= A:sin6'cos6' — Ml COS0 + U2 sin^ (11) 
— cot 0(ui sin (/) + U2 cos (/)). 
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It can be simplified by using the controls vi and V2 such that: 

Vl = -wi COS0 + U2 sin0 .^^ 
V2 = —ui sin — U2 cos (p. 

Note that the cost C is not modified since vl + = u\ + u^. The Hamiltonian 
H reads now: 

H = — fcrcos^ Bpr + pe{ksa\9 cosB + vi) ~ V2P4, cot9 — 2(^1 + '^2)- (13) 
Using the maximization condition, one arrives at: 

Vl ^ pg, V2 ^ -cot 9p4,. (14) 
Plugging the expressions of vi and V2 into H leads to: 

H = -fcr cos^ 9pr + kcos9sin9pg + ^pl + i cot^ 9p'^. (15) 

The Hamiltonian dynamics is then characterized by the following differential 
equations: 

= —krcos^9 
— k sin 9 cos 9 + pg 



cot^ 



Pr = kpr COS 9 

pg = -krprSm{29)-kcos{29)pe+pl^ 

V6 = 0. 



(16) 



In the STIRAP process, the coordinate X2{fy remains equal to for any time 
i, which avoids the dissipation effects and maximizes the transfer to the state 
3). In the spherical coordinates, we have 9(fy = 7r/2 and 9(t) ~ which leads 
to pg(i) = (/'(t) = 0, i.e. no motion is possible. This remark is an illustration of 
the well-known fact that the exactly adiabatic trajectory is singular and needs 
an infinite transfer time. 

2.3 Geometric analysis of the Hamiltonian dynamics 

This preliminary analysis does not give any insight into the geometric character 
of the adiabatic solution, and a more intricate understanding of the Hamiltonian 
dynamics is required. This can be done first by proving that the Hamiltonian 
given by Eq. is Liouville integrable [38J, i.e. it has as many constants 

of the motion as degrees of freedom. Since H does not depend on we can 
check that p^ is a constant of motion. Another constant of the motion can be 
found by noting from Eq. (jl6l) that the product rp^ is constant. Using the 
change of coordinates p = logr and introducing the new set {p^pp) of canonical 
coordinates [33, one deduces that the Hamiltonian iJ becomes: 

B = - kp pcos^ 9 + kain 9 cos 9pg + ^{pj + cot^Opl), (17) 

and does not depend on p. The functions H, p^ and pp are three constants of 
motion and the system is therefore Liouville-integrable. 



5 



To simplify the description of the Hamiltonian dynamics, we will consider in 
the following only a two-dimensional degree of freedom system in the coordinates 

(/)) taking pp as a parameter. This is possible since the dynamics does not 
depend on p. The Liouville-Arnold theorem tells us that a given trajectory 
defined by some initial conditions will curl around a torus in the phase space 
|381I39| . This torus can be either regular or singular, the former structure giving 
periodic solutions while the latter leads to trajectories with an infinite period. 
For one-dimensional Hamiltonian systems, note that the singular tori are the 
well-known separatrices. Hence, one could expect that the STIRAP solution is 
linked to such singular tori, which is therefore crucial to characterize. For that 
purpose, we consider the so-called energy-momentum map: 

i.e. the set of all the possible values of the constants of the motion H and 
P0, which is parameterized by the constant pp. A rigorous application of the 
Liouville-Arnold theorem gives that the inverse image J-'~^{H,pcj,) of the point 
{H,p^) of the energy-momentum set is a two-dimensional regular torus if the 
two gradient vectors of the phase space V-ff = ( ^| , , ^ , ^ ) and Vp^ = 
(0,1,0,0) are not parallel for any point of the torus. If the two vectors are 
parallel for some points, then the pre-image is no longer a regular torus, but 
a singular one. The topology of the singular torus can be of different types: a 
point (for an equilibrium), a circle (for a periodic orbit) or a more complicated 
two-dimensional structure. This latter geometry can be determined from the 
singular reduction theory, which is applied for this example in the appendix 
El Here, we only compute the position of the different singular points in the 
energy-momentum diagram. If the two gradients ViJ and Vp^ are parallels 
then the matrix defined by 

/ k sin 9 cos + pg 

cot^ 9 Pel, 
kpp sin(2^^) + k cos{2e)pe - 
V 

is not of full rank, i.e. equal to 2 [39]. We obtain two different cases. In the 
first one, which corresponds to the STIRAP process, we have {0 = 7r/2,pe = 0) 
and H = 0. In the second situation, one arrives at: 

P0 = —ksmOcosO 

P^ = ±^'^{kpp sm{20) + k cosi20)pg), 

which gives the equation of the boundary of the energy- momentum set. The cor- 
responding energy-momentum diagram is represented in Fig. [TJ Its qualitative 
structure does not depend on the value of pp. 

As mentioned above, a singular reduction |52] allows to make more precise 
the nature of this singularity (see the appendix |^ for details) . In the energy- 
momentum diagram depicted in Fig. [U each point of the horizontal singular 
line of equation _ff = is associated to a bitorus in the phase space, i.e. two tori 
glued together along a singular circle. Straightforward computations show that 
the points of this circle satisfy X2 = Px2 = Oj i-e- the equation of the STIRAP 
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(18) 
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Figure 1: (Color Online) Image of the energy- momentum map J- (gray) for the 
Haniiltonian of Eq. (jlSp . The singular points are represented by solid blue 
(black) and red (dark gray) lines. The red line indicates the positions of the 
image of the bitori (see the text for details). Two trajectories are schematically 
represented on a regular and on a singular tori (the dots represent the initial and 
final points of the trajectory). The upper trajectory represents a standard oscil- 
lating solution obtained with the energy-minimum cost. The lower trajectory is 
the ideal adiabatic trajectory, that can not be reached with the energy-minimum 
cost. Numerical values are taken toh^ — k — \. The different quantities are 
unitless in all the figures. 
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Figure 2: (Color Online) Top: Projection in the {x2,Px2) plane of a trajectory 
on the bitorus for iJ = 0, = 0.1, Pp = 1, in the case where 6(0) ^ 0. The 
parameter k is fixed to 1. Bottom: Projections onto the sphere of different 
trajectories lying on a bitorus. The green (light gray) and red dots (dark gray) 
are respectively the initial and final positions of the trajectories. The great circle 
represents the quantum states with no population in the intermediate level |2). 
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Table 1: Characteristics of different extremal solutions. The Freq. and Amp. 
columns indicate respectively a frequency and an amplitude averages of the 
fields ui and U2. 



H 


P0 


T 


Freq. 


Amp. 


C 


^1{T) 





O.f 


4 





1.7 


2.7 


10-7 


0.33 


15 


30 


1.5 


1.26 


33.9 


0.82 


0.4 


45 


8 


1.5 


4.23 


73.6 


0.93 


4.6 


30 


10 


50 


9.81 


94.3 


0.98 



solution. For comparison, two schematic trajectories on a regular torus and on 
a bitorus have been plotted in Fig. [1] We notice the difference of structure 
between these two solutions. The set of all the trajectories belonging to a given 
bitorus can also be characterized. Two different cases can be distinguished. In 
the first one, = -j, then pg = Q and the system is on the singular circle which 
means that no motion is possible. In the second subset, 6{t) ^ at some time 
t, but we observe numerically that the trajectory returns at longer times to the 
singular circle such that 6* = f- This point is illustrated in Fig. [2l 

The last step of this preliminary study consists in solving numerically the 
optimal control problem. For that purpose, we choose some initial values of the 
momenta pp, pg andp^ and we propagate the Hamiltonian equations (jl6p during 
the time T. The distance to the target state is computed at this final time. An 
example of solution is represented in Fig. [31 Note the oscillatory behavior of 
the extremal solution, which is far from the expected smooth and monotonic 
evolution of the adiabatic pulse. These oscillations do not allow to limit the 
losses from the intermediate state and the final population of the state |3) is 
only of 0.82. Such a solution belongs to a regular torus close to the singular 
one. In spite of this proximity in the energy-momentum diagram, two completely 
different behaviors are obtained. In addition, a systematic numerical analysis 
shows that the frequency of the oscillations increases if the control amplitude 
increases as can be seen in Tab. [TJ We observe that the population transfer is 
better, but the control frequency is larger. On the other hand, if one tries to 
reduce the oscillatory character of the solution by staying on the singular line 
or close to it then it is not possible to reach the target state with efficiency, as 
illustrated by the first line of Tab. [TJ We thus conclude that the limit of high 
energy does not allow to recover the STIRAP process. The other intuitive limit 
would be a longer control duration, but this is not favorable to STIRAP, as 
can be seen by comparing the second and third lines of Tab. [1] The STIRAP 
process is therefore not an intrinsic optimal solution of this problem. A peculiar 
cost functional has to be chosen in order to highlight the optimal properties of 
the STIRAP technique. 

2.4 The STIRAP cost 

The study presented above makes clear that the energy minimum cost is not 
well suited to the STIRAP process. The objective of this section will be to 
determine a cost, which could force the optimal solution to follow the STIRAP 
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Figure 3: (Color Online) Projection onto the sphere of an optimal trajectory for 
the energy minimization control problem (top panel) and of the corresponding 
fields (bottom panel) . The pump and the Stokes fields are respectively depicted 
in blue (black) and in red (dark gray). The parameters are fc = 1 and T = 30. 
The values of the adjoint states are = 15, Pp = 69 and the Hamiltonian H 
of Eq. ^ is equal to 0.33. 
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trajectory. From a dynamical point of view, a basic argument is to prevent the 
Hamiltonian trajectory to curl around a torus. This goal can be achieved by 
reducing the oscillations around X2 = 0, which leads in spherical coordinates to 
the cost: 



C = 



[k sin 9 cos 9 + viYdt. 



(19) 



This choice also avoids populating the state |2) as shown below. The Hamilto- 
nian of the PMP has the form: 



H = —kr cos 



k cos 9 sin 9p0 + vipe 



—V2 cotOp^ — -{ksin9cos9 + wi)^. 

The cost C depending on vi and not on W2, we choose to optimize the Hamil- 
tonian H only with respect to vi . Optimizing also V2 would lead once again to 
a static solution. We obtain vi — pg ~ fcsin^cosfl and plugging this expression 
into the Hamiltonian H gives: 



H = —krpr cos^ 9 -I- -pg — V2 cot ( 



(20) 



where V2 is a given time-dependent function. The equations of motion become: 
f — ~krcos^9 



(j) = — cot 9 V2 

Pr = kpr COS^ 9 

pe = —2krprSm9cosl 

Vd, = 0. 



(21) 



p<f, 



^V2 



From Eq. (|2ip , it is clear that the solution pg (t) = for any time t minimizes 
the cost C. As a consequence, we get that 9{t) is constant along this trajectory 
and the expression of the second control field V2 comes from pg — 0: 



V2 



2krpr sin 9 cos 9 



(22) 



Here the functions rpr, p^ and 9 are constants of motion, V2 is also constant. 
Using Eq. ([22)). one gets the extremum value of H, which reads as: 



H = 2krpr cos^ 9{sm^ 



(23) 



The STIRAP solution can be recovered from the limit case -> 7r/2. In this 
limit, we get H ~ Q and vi ~ 0, leading to ui coscj) — U2 sincj) and 



U2 _ Xi 
Ui X3 '' 



(24) 



which is characteristic of a STIRAP sequence. 

To numerically solve this control problem, one should also note that there 
exist some constraints on the initial momenta. Indeed, if we denote by t the 
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Figure 4: (Color online) Optimal control fields obtained with the Stirap cost 
(bottom) and the projection of the corresponding trajectory onto the sphere 
§ = 1 (top). In the bottom panel, the Stokes and the Pmnp pulses 
are respectively depicted in blue (black) and red (gray). Numerical parameters 
are taken to be fc = 1 and T = 80. The initial momenta are pg = 0, = 0.1, 
Pp = 10. The Hamiltonian H of Eq. ^ is equal to 10"^. 
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typical relaxation time, the following condition holds on the final time: T << t. 



As T = 



and T = 



2v2 cot 6 



this leads to the relation: 



< 1. 



Arpr sin f 

An example of optimal trajectory is given in Fig. U) Note the counterin- 
tuitive order of the two control fields, one of the main features of a STIRAP 
process. We can therefore conclude that a STIRAP-like solution can be designed 
from the PMP if the right cost is considered. 



3 Extension to a four-level quantum system 

We apply the same reasoning to a four-level quantum system with a tripod 
structure [371[ini[lIllllll3[71[Ml[M]as illustrated in Fig. [3 




Figure 5: (Color Online) The four-level tripod quantum system, ui, U2 and uj, 
stand for the pump and the two Stokes pulses respectively, k is the dissipation 
rate of the second level |2). 



The dynamical evolution of this system is described by the Hamiltonian: 

/ -ui \ 



H 



(25) 



ui —ik —U2 — U3 
M2 

\ U3 y 

where ui, U2 and M3 stand for the pump and the two Stokes pulses respectively. 
The aim of the control is to transfer the population of state |1) to a superposition 
of states |3) and |4), while preventing the system to populate the state |2). We 
first consider the energy minimization cost C given by 



[ul{t) +ul{t) +ulit)]dt. 



The pseudo-Hamiltonian of the PMP reads: 

H = -~kx2Px2 + ui{xipx. 

+U3.{X2PX4 - X4.PX2] 

Introducing the spherical coordinates: 



X2PxJ - 



U2{x2Px^ ~ X3PX2) 

-ul{t)+ul{t)]. 



(26) 



(27) 



: r cos 9i sin 62 
X2 = r cos 02 
x^ = r sin 9i sin 62 cos 83 
X4 = r sin 01 sin 62 sin ^3 
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the dynamical system takes the form: 
f = — fcrcos^02 

91 = ui sm cot 62 + U2 '^^ ^ 
+U3 cot 02 cos 01 sin 03 

92 — k sin 02 cos 02 — ui cos 0i + U2 sin 0i cos 03 ' 
+U3 sin 01 sin 03 

g cos 6*2 sin 6*3 , cos 62 cos 83 

73 — "2 sine^ sinea 3 sinfli sinfla 

Using the foUowing rotations on the control fields: 

V2 = U2 COS 03 + Us sin 03 
W3 = — 7i2 sin 03 + U3 COS t 



(29) 
(30) 

i7 = — r cos^ 02Pr + COS02 sin02P02 + ''^i cot 02Pei 



3 

and 

wi — Til sin 01 + t>2 cos 01 
u'2 — —Ml cos 01 + U2 sin 01 

the Hamiltonian H becomes: 



(31) 



The analysis of the different Hamiltonian trajectories used in Sec. [2]can be done 
along the same lines. It can be shown that the Hamiltonian H is integrable 
since it possesses four constants of the motion: H, Li = x^P^^ — x^Px^, L3 — 
xiPx^—XiPx^ and L4 = xiPx^ —x^Px^ ■ Here again, we obtain that the adiabatic 
trajectory lies on a singular torus in the phase space. In this example, we have 
to deal with a four dimensional torus in a eight dimensional phase space, which 
prevents any three dimensional picture. As in the STIRAP case, the optimal 
solution presents an oscillatory structure in the 02 variable which is not present 
in the adiabatic solution. Following the same intuition as in Sec. ^ we give up 
the minimization of the energy and we choose a cost that allows to minimize 
the oscillations around 02 = f: 

C = / e^dt^ I {ksm02Cose2 + W2)^dt. (32) 
Jo Jo 

Plugging this cost into the Hamiltonian and optimizing with respect to 'W2, 
which is the only control in the cost, we get: 

W2 — Pe2 " k sin 02 cos 02 , 

and the Hamiltonian finally reads: 

H = -krpr cos^ 02 + ^pj^ + wi cot 02Pei 

cotSape, (33) 

It is then straightforward to check from the equations of motion that the solution 
Pe2 (t) — Oj for any time t, minimizes the cost C. One deduces that 02 is constant 
along this trajectory and we obtain a relation between V3 and wi from pe^ ~ 0: 

V3 — {2krpr cos 02 sin'^ 02 — wipoi ) ■ (34) 
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Figure 6: (Color Online) Example of an optimal solution mimicking the adia- 
batic evolution and leading to a superposition of the states |3) and |4). The 
top panel depicts the evolution of the population of the different states |1), |2), 
|3) and |4) respectively in blue (dark gray), red (light gray), black and green 
(dashed line). In the bottom panel, the Pump, Stokes 1 and Stokes 2 pulses 
are represented in blue (dark gray), red (light gray) and black, respectively. 
Numerical values are pp = 100, wi = fc — 1. 
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Hence, by tuning the values of wi, which is assumed here to be constant, and 
of the initial momenta, we can reach any final superposition of the states |3) 
and |4). For example, by taking wi — 1, rpr — 100, pei{0) — 16.85, psaCO) = 
0,^^3(0) — —1, we obtain the same final population in the states |3) and |4) as 
shown in Fig. [HI Note the counterintuitive order between the pump and the 
Stokes pulses. The corresponding adiabatic process is associated to the limit 
O2 ^ 7r/2. 

4 Conclusion and open questions 

We have reached the main goal of this study, which was to exhibit a connection 
between optimal control theory and adiabatic techniques in quantum systems. 
A geometric analysis of the Hamiltonian dynamics constructed from the PMP 
gives a complete overview of all the optimal solutions for a given cost. For 
the energy minimization problem, we have shown that the adiabatic pulse can 
be associated to a peculiar Hamiltonian singularity of the problem. The adia- 
batic solution being a singular Hamiltonian trajectory, it cannot be approached 
smoothly by any optimal control field. Therefore, an adapted cost functional 
has to be chosen to enforce the optimal solution to follow the structure of the 
adiabatic pulse sequence. This study done for three and four level quantum sys- 
tems is expected to be generalizable to more complex quantum dynamics where 
an adiabatic control scheme can be used. On the theoretical side, this study is 
also the first step in the understanding of the role of Hamiltonian singularities in 
control processes. Although much more work need to be done to advance in this 
entirely new field, it seems a promising way in order to design optimal control 
fields with robustness properties. For instance, it would be interesting to apply 
the same idea to other standard adiabatic processes such as the frequency chirp 
in a two-level quantum system. 
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A Singular reduction theory 

In this appendix, we apply the singular reduction theory |39| in order to deter- 
mine the nature of the singular tori presented in Sec. [5] We refer the reader 
to Ref. [445 for a pedagogical introduction to these tools. Singular reduction 
theory is a general technique that gives a global overview of all the Hamilto- 
nian trajectories and of the corresponding geometrical structures in which they 
evolve. This theory is roughly based on the reduction of the dimension of the 
phase space by making use of a constant of the motion. However, caution should 
be exercised due to the existence of singularities in the problem. The applica- 
tion of this general procedure to the example of this work can be summarized 
as follows. 

To simplify the analysis, we work on the dynamics projected on the sphere, 
the constant r -pr playing the role of a parameter. In cartesian coordinates, this 
leads to the following constraints: 

xl + xl + xl = l and xip^^ + X2PX2 + xsPxa = rpr. (35) 
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We consider the flow of the constant of motion = xiPx^ — x^px^ to reduce 
the dimension of the problem. We first introduce the six invariant polynomials 
under the flow of p^, which form a basis for the polynomial functions V which 
Poisson-commute with p^, i.e. V = {V,Pcf,} = 0. It is straightforward to check 
that the following polynomials are invariant |39j: 

TTl = X2 
71'2 = Px2 

71-3 = XlPx, - X3Pxi (= P,p) .gg., 
■71^^x1+ xl 

TTe = xipx^ + xzpx3 

These polynomials obey by construction to the reduced phase space equation: 

TTg + TTg = 7r47r5. (37) 

Since the Hamiltonian H defined in Eq. (jl5D commutes with p^, it can be 
expressed as: 

H = -A:7ri7r2 + - {-kIt^a + TTaTTs - 2nin2ne) (38) 
Using the two constraints of Equations (P7)) and ([M)) become: 

(rpr - 7ri7r2)^ + TTg = 7r4(l - TTj), (39) 

and 

H = -{k + rpr)TTin2 + - {TrfiTA + nlnf + tt^) . (40) 

For a given value of 7r3(= p^), we get two surfaces in the space (tti, 7r2, 714). The 
dynamics takes place at their intersection. An example of such intersection is 
plotted in Fig. [21 This eight-like shape is the mark of a specific type of singular 
tori, the bitorus which is represented in Fig. [31 
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